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Abstract 

With sufficiently liglit up and down quarks the isovector (oq) and isosinglet (/o) scalar meson 
D , propagators are dominated at large distance by two-meson states. In the staggered fermion for- 



mulation of lattice quantum chromo dynamics, taste-symmetry breaking causes a proliferation of 



^ ■ two- meson states that further complicates the analysis of these channels. Many of them are unphys- 

, ical artifacts of the lattice approximation. They are expected to disappear in the continuum limit. 

CN \ The staggered-fermion fourth-root procedure has its purported counterpart in rooted staggered 

(•~^ ■ chiral perturbation theory (rSxPT). Fortunately, the rooted theory provides a strict framework 

. that permits the analysis of scalar meson correlators in terms of only a small number of low energy 

. ^ couplings. Thus the analysis of the point-to-point scalar meson correlators in this context gives a 

X' 

' useful consistency check of the fourth-root procedure and its proposed chiral realization. Through 

numerical simulation we have measured correlators for both the oq and /q channels in the "Asq- 
tad" improved staggered fermion formulation in a lattice ensemble with lattice spacing a = 0.12 
fm. We analyze those correlators in the context of rSxPT and obtain values of the low energy 
chiral couplings that are reasonably consistent with previous determinations. 

PACS numbers: ll.15.Ha, 12.38.Gc, 12.39.Fe, 14.40. Cs 
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I. INTRODUCTION 



The recent evident successes of numerical simulations of QCD with improved staggered 
fermions demand a thorough examination of its most controversial ingredient, namely, using 
fractional powers of the determinant to simulate the correct number of quark species (the 
"fourth root trick"). The procedure is known to introduce nonlocalities and violations of 
unitarity at nonzero lattice spacing [1]. If these problems do not vanish in the continuum 
limit, they may even place the theory in an unphysical universality class. There are, however, 
strong theoretical arguments [2, 3, 4, 5] that the fourth-root trick is valid, i.e. that it produces 
QCD in the continuum limit. 

One may also test the fourth-root procedure numerically. One can, for example, check 
that taste symmetry gets restored as the lattice spacing gets smaller, by looking at the 
eigenvalue spectrum [6, 7, 8, 9, 10], the Dirac operator [11], or the pion spectrum [12]. Al- 
ternatively, low-energy results of staggered fermion QCD simulations can be compared with 
predictions of rooted staggered chiral perturbation theory (rS^PT) [13, 14]. Since staggered 
chiral perturbation theory becomes standard chiral perturbation theory in the continuum 
limit, agreement between rooted QCD and (rS^PT) at nonzero lattice spacing would sug- 
gest that, at least for low energy or long-range phenomena, lattice artifacts produced by the 
fourth root approximation are as harmless as those produced by partial quenching. Partial 
quenching also induces unitarity violations, but they disappear in the limit of equal valence 
and sea quark masses. 

There are two recent tests of agreement between rooted staggered fermion QCD and 
rSxPT: (1) Measurements of the light pseudoscalar meson masses and decay constants in 
partially quenched and full staggered fermion QCD fit well to expressions derived from 
rSxPT [15]. A byproduct of this fit is a determination of the low energy couplings of the 
theory. (2) The topological susceptibility measured in full QCD agrees reasonably well with 
predictions of rS^PT [16]. 

In the present work we examine scalar meson correlators in full QCD and compare their 
two-meson content with predictions of rS^PT. Since the appearance of the two-meson 
intermediate state is a consequence of the fermion determinant, an analysis of this correlator 
provides a direct test of the fourth root recipe. The channel has been studied recently 
in staggered fermion QCD by the MILC collaboration [17] and UKQCD collaboration [18]. 
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Both groups found that the correlator appeared to contain states with energies well below 
possible combinations of physical mesons. 

A simple explanation of the nonstandard features of the scalar correlators is provied by 
rSxPT [19, 20]. In that theory all pseudoscalar mesons come in multiplets of 16 tastes. The 
pattern of mass splittings is predicted by the theory. The tt and K multiplets are split in 
similar ways. The rj and rj' mesons, however, are peculiar, because their masses are shifted 
by the axial U{1) anomaly. Since the anomaly is a taste singlet, only the taste singlet r] 
and 7]' acquire approximately physical masses. Some of the remaining members of the r] 
multiplet remain degenerate with the pions. According to taste symmetry selection rules, 
any two mesons coupling to a taste-singlet Oq must have the same taste. But all tastes 
are equally allowed. Among other states, the taste singlet couples to the Goldstone 
pion (pseudoscalar taste) and an 77, also with pseudoscalar taste and of the same mass. 
This spurious two-body state at twice the mass of the Goldstone boson accounts for the 
anomalous low-energy component in that channel. 

This explanation raises concerns. Clearly, only the taste singlet 77 approximates the 
physical state, since it is the only member of the multiplet subject to the anomaly. So if 
the other ?7's are not allowed as external states, we have violated unitarity in the sense that 
some intermediate states are not allowed as external states. Further examination of the 
taste multiplets in the intermediate states reveals that in addition to the several unphysical 
Tir] taste combinations, there is a negative norm "ghost" contribution in the taste singlet 77 
meson leg [21]. Remarkably, all lattice artifacts resolve themselves in the continuum limit, 
however. The taste multiplets become degenerate, the two-body states merge, and the 
ghost state cancels the spurious taste combinations, leaving only the taste-singlet mesons. 
To achieve this cancellation requires following the rules of flavor counting in rS^PT. 

In the present work we extend the analysis of Ref. [19, 20] and carry out a quantitative 
comparison of measured correlators and predictions of rS^PT. Despite the considerable 
complexity of channels with dozens of spectral components, the chiral theory models the 
correlators precisely in terms of only a small number of low energy couplings, which we may 
determine through fits to the data. 

This article is organized as follows. Following a review of some needed results from S^PT 
in Sec II, we derive the chiral predictions for the ao and /o in Sec III. We present results of 
our fits to the predicted forms in Sec IV and conclude in Sec V. 
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II. ELEMENTS OF STAGGERED CHIRAL PERTURBATION THEORY 



In this section we give a brief review of rooted staggered chiral perturbation theory with 
particular emphasis on the tree-level pseudoscalar mass spectrum. We obtain the rooted 
version of the theory through the replica trick, according to which each quark flavor, u, d, 
and s, comes in four tastes and is repeated times [22]. We calculate various quantities 
in the replicated theory, and in the final step, we set rir = 1/4 to obtain the correct flavor 
counting. 

The low energy effective chiral theory is formulated in terms of the meson field 

16 

* = E T^Ty (1) 
b=i ^ 



where = {1, Cs, "^^sC/^ • • • } are Dirac gamma matrices and </>'' is a Sn^ x Sn^ matrix with 
rows and columns labeled by the flavor and replica index ur, dr, and sr. The staggered 
chiral action is written in terms of the unitary matrix S = exp(2z$//): 

5(S,m)= J d'y i^^TT{d,J:^d^^^) - i^TiiM^^ + M^J:) + ^4>l, + a'ViJ:)^ . (2) 

The low energy couplings at this order are /, /i, and the quark mass matrix Ai = It 
/rdiag(mu, m^, mg), where It is the unit matrix in taste space and is the unit matrix in 
replica space. The axial anomaly appears through the mass term ttIq. It involves the flavor- 
singlet taste-singlet field 00/ = J2f,r 4>fr,fr/ V^^-- The taste-breaking term V is a linear 
combination of operators [13, 14, 23] 



(3) 



where 



O3 

0,v 



TrlTo.STo.St 



[Tr(To,^S)Tr(To,^S) + h.c] 
[Tr(To,^5S)Tr(To,5^S)+h.c.: 



4 
1 

4 

^[Tr(To,^STo,^S)+h.c. 



[Tr(To,^5STo,5/,S) + h.c. 



Tr(To,^E)Tr(To,^S^ 



(4) 
(5) 
(6) 
(7) 
(8) 
(9) 



OsA = \ [Tr(ro,^5S)Tr(To,5^St)] (10) 
a = E Tr (To,^.STo,,^St) . (11) 



Without the anomaly and taste-breaking term the tree-level masses of the pseudoscalar 
mesons with quark flavor content /, /' cirG , clS usual, 

M]j,^, = fi{mf + mf,). (12) 

The taste-breaking term splits the nonisosinglet states (vr;, and Kb) to give 

where to leading order the multiplets split five ways, 



As = 
16 

A^, = p(2C3 + 2C4 + 4C6) 



A^5 = 7i(Ci + 3C3 + C, + 3Ce) (14) 



A^ = ^(Ci + C3 + 3C4 + 3C6) 



A, = ^(4C3 + 4C4), 

which we label P, A, T, V, and J, respectively. This predicted multiplet pattern has been 
well confirmed in simulations [17, 24]. 

We will be working with degenerate u and d quarks (m„ = rrid = m^), so it will be 
convenient to introduce the notation 

M^fe = 2/im^ + a^Ab 

= 2fims + a^Ab (15) 
= K^^i + ^s) + a^Ab 

The isosinglet states [rj and rj') are modified both by the taste-singlet anomaly and by 
the two-trace (quark-line hairpin) taste- vector and taste-axial- vector operators 02V-, O2A, 
O5V and O^A- When ml is large, in the taste-singlet sector we obtain the usual result 

Ml = Im'^j + ImI, 

M,,j = 0{ml). (16) 



In the taste- axial- vector sector we have 

M^,A = \\MIa + ML + ^UrdA - Za] 

Ml A = \ [Ml A + + ?,n,.5A + Za] (17) 
Zl = {MIa - Mlj^f - 2nr5A{Mlj, - M^^) + <dnl5% 

where 5a = a^S'^^ = a^l6{C2A — C5A)/ P, and hkewise for A ^ V. 

In the taste-pseudoscalar and taste-tensor sectors, in which is there is no mixing of the 
isosinglet states, the r], and r]', by definition have quark content (uu + dd)/^/2 and ss, 
respectively, and masses 

Ml, = Ml, ■ Ml, = Ml, (18) 

In Table I we list the masses of the resulting taste multiplets for the lattice ensemble used 
in the present study with taste-breaking parameters Sa and 6v determined in Ref. [15, 17]. 

b TT, K, T], rl, 
P 0.1594 0.3652 0.1594 0.4927 
A 0.2342 0.4036 0.1843 0.5129 
T 0.2694 0.4250 0.2694 0.5384 
V 0.2966 0.4428 0.2825 0.5491 
1 0.3205 0.4592 0.4958 - 

TABLE 1: Masses of pseudoscalar meson taste multiplets in lattice units for the MILC coarse 
(a = 0.12 fm) lattice ensemble /? = 6.76, arriud = 0.005, omg = 0.05, as measured or inferred from 
measured masses and splittings. The mass of the 77^ depends on the anomaly parameter ttt-q. 

III. SCALAR CORRELATORS FROM S^PT 

In this section we rederive the "bubble" contribution to the channel of Ref. [20], using 
the language of the replica trick [13, 25], and then extend the result to the /o channel. 

We match the point-to-point scalar correlators in chiral low energy effective theory and 
staggered fermion QCD by matching the Green's functions, which are defined through the 
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generating functionals of the respective theories: 

92 \ogZ 



(19) 



dmfj>{y)dme',e{0) 

For this purpose the quark mass term diag(mu, ma, mg) is converted to a local meson source 
rriffily) (including flavor off-diagonal terms) in both S^PT and QCD. 



A. Scalar correlator in staggered fermion QCD 

First we review the construction of the needed correlators in staggered lattice QCD, where 
the generating functional is 

Z{mff,)= J dUexp[-Sg{U)]det[M{U,mff)Y/^ . (20) 

Here U are the gauge link variables, Sg{U) is the gauge action, and M is the fermion matrix 
including flavor components. We work on a lattice of spacing a and dimension x Nt and 
label sites by the integer four- vector x^. Hypercubes of size 2^ are similarly labeled by y^, 
so = 2?/^ + ?7^. 

Staggered fermion meson correlators can be defined in the one-component basis of the 
Grassman color vector field Xf{^) or in the spin-taste basis of the field qj'^iy) with spin label 
a and taste label a. The fields are related through 

o 

Xf{2y + V) = 2Tr[rtg^(y)] (21) 

where = 70" 72^73^? and the sum over r] runs over sites in the 2^ hypercube labeled by 
y. The lattice y has spacing A = 2a. 

For constructing the meson correlators via the functional derivative (19) we need to 
introduce the source term into Lagrangian 

Sm = a^Y.Xf{x)xf'{x)mfj>{x) . (22) 

X 

To express the source in term of the spin-taste basis we use the relation 

a%i2y + v)Xf'i2y + r/) = ^ C(r, v)Pfj'Ay) (23) 
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with 



and 



C(r,r/) = Tr(rtrtr,r)/4 (24) 



p/j',r(y) = g/(y)r®rg^,(y) . (25) 

The direct product Ts ® acts on spin and taste components, respectively. So we obtain 

Srn = A'^Yl Pfj'Ay)i^fj'Ay) (26) 

y,r 

with 

rrifj'Ay) = ^ E C(r, v)mfA2y + r/) (27) 

The desired source for the scalar density, mfj/j{y), has T = I and C{I,ri) = 1. It is the 
component of mfji{x) that is constant over a 2^ hypercube. The other terms mjj'_r(l/) are 
sources for the other local staggered mesons. 

A particular correlator is obtained by differentiating the generating functional with re- 
spect to the appropriate source mass terms. The general two-point function is, then, 

92 logZ 



(9m/j/,r(?/)ame',e,r'(0) 



= A«(p;j,r(y)Pe',e,r'(0)). (28) 

m^ji (x)=5 J jiruf 



The above quantity will be calculated for T = I also within S^PT below. 

Now, we need to relate the quantity (28) to the correlator generated from the code. In 
practice the simulated correlator is computed from a point source at the origin 

Oe,e',src = o'Xe' (0)Xe (0) = ^ EPe,e'r(0), (29) 

8 r 

and a single time-slice sink operator at time t = 2t + tiq, 
0/,/'sink(y,r) = a3^x/'(2y+r/,r)x/(2y+r/,r) = A'[pfj,j{2y,t) + {-rpfj,,o5i2y,t)], (30) 

V 

where we have used relation (23), Ci^^v) = 1 a-^d ({05,7]) = {—)^°- Note that the sink 
operator is defined on spatial cubes y but all time slices r. 
In this language the computed correlator is 

Cf,f';e,e'{P, ra) = J2 exp{ip- yA) (O fj'^i^y,{y, r)Oe,e'src) • (31) 
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The meson taste is conserved, so the correlator separates into nonoscillating and oscillat- 
ing components for a taste-singlet scalar contribution and a taste-axial-vector pseudoscalar 
meson contribution, respectively. 

Cfj'.e,e'{p, to) = Cfj'.e,e',l{p, TO) + (-)^C/j/;e,e',05(P, ™), (32) 

where 

Cfj'-e,e',T{p, ™) = X] exp{ip- yA) {pfjf^rC^y, i)Pe,e',r(0)) . (33) 

y 

The correlator has a quark-line-connected part and may also have a quark-line discon- 
nected part: 

C /J';e,e',r(P, ra) = C /,/';e,e',r,conn(p, ™) + C /,//;e,e',r,disc(p, ™) (34) 

The quark-line disconnected part appears only in the taste-singlet isosinglet correlator. 

To compute the correlator we need to express it in terms of quark propagators. So we start 
from the definition of the correlator in Eq. (31), substitute the definitions of the operators 
in Eqs. (29) and (30) and use the relation 

52 log Z 



a {xfC^y + v)xf'{'2y + v)xe{v')xe'{v')) 



dmfj>{2y + r])dme',e{v') 



(35) 

m.j,f,{x)=Sff,mf 



^ E C(r, r/)C(r', v') {Pfj'Ay)Pe',eAo)) 



which follows from identity 



d 1 d 

dm{2y + r]) " 16 ^ dmr{y) ' ^^^^ 



Finally we arrive at the point-to-point correlators 



C/,/';e,e',conn(p, ™) = - exp(zp ■ xa) {Tt[Mj\x, t; 0, 0)My^\x, r; 0, 0)]) 5e/5e'/' 

C/,/';e,e',disc(p, ™) = ^^exp(zp- xa) {^Tr [M^^ (f , T ; f , r)]Tr [M^"^^ (0, 0; 0, 0)] ^ (5ee'5//' , 

X 

where we have used Eq. (20) and the normalization M = 2D + 2am for the Dirac matrix. 
We keep the momentum p small, so we can neglect variation of the exponential over the 
hypercube. 

As it is computed, at zero momentum the quark-line disconnected correlator includes the 
vacuum disconnected piece: 

Cj,e,o = ^ (Tr[M7i(0, 0; 0, 0)]) (Tr[M-^(0, 0; 0, 0)]) . (37) 
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B. Scalar correlator in S^PT 



The continuum generating functional for scalar correlators in S^PT 

^sxpt("1//') = J [dS] exp[-5'(S, m///)], (38) 

where S{'E,mffi) is given by Eq. (2). We do not include explicit scalar meson fields in the 
chiral Lagrangian, but add their contributions in the final expressions. To match the func- 
tional derivatives (19) we approximate the continuum integration in the chiral theory with 
a sum over hypercubic volumes of size and differentiate with respect to a constant source 
inside that volume. In this case mfjijiy) = mf f/{y). The source is also constant over repli- 
cas of the same flavor. The space-time volume equals that of QCD, namely, A'^{L/2)^Nt/2 
for A = 2a. We use the integer four vector y to label the hypercubes in the chiral theory. 
The functional derivative in S^PT is 

dmfj'{y)dme',e{0) 7t>^ hrj'r H > ) er' ,e'r' / ^ ' 

At tree level the action (2) has no explicit quark- ant iquark scalar meson fields, but it gener- 
ates the two-pseudoscalar-meson "bubble" terms in the correlator. The functional derivative 
(39) corresponds to (28) with F = J. We use B to denote the bubble contribution corre- 
sponding to the correlator (33) 

BfJ'-e,e',l{P,tA) = — ^exp(ip-yA) (p/,j/j(?/)pe,e',/(0)) (40) 

y 

y ^i*"' 
We introduce its time Fourier transform 

Nt/2 

Bfj>.e,e',i{p) = ^ exp{ipotA)Bfj>.e,e'j{p,tA). (41) 

t=0 

At tree level the vacuum expectation value reduces through Wick contractions to products 
of meson two-point functions. In momentum space we have, generically, the Euclidean 
correlator 

^'^^y^'^^^^^ = A^{L/mN^/2) ? '^P^'^ ■ (42) 
where {(f){-k)(j){k)) = l/{e + m^). So 

/A( \Afn\\ 1 ^ exp[-E{k)tA + ik-yA] 

(<Piy)m) = E (43) 
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for E{k) = y\k\'^ + m'^- In terms of momentum components, the general term in the corre- 
lator becomes 

Bf,f';e,e'AP) = -J-T.^MW-yA) J2 E (</'/r,^;s(^.^)</'5s,/v(^,^)C',,'.'(0)<'.Vr'(0) 

y g,s,r,b g' ,s' ,r' ,b' 

= 8(L/2)3(iV,/2)A2?^£^,5,, Mr,si-k)€',Mk)) 

We have used the fact that the bubble term, by definition, does not include the vacuum 
disconnected piece corresponding to Eq. (37). 

There are two types of two-point functions, namely, the connected two-point function for 
all tastes: 

and the additional disconnected contribution for the taste-singlet, taste-axial-vector, and 
taste-vector mesons. For the taste singlet it is 



"I^Ufri k)(l)f,r',g's'ik)) - '" J'' ' ^ ^ . (46) 



Here we have already decoupled the taste-singlet rj' by taking ml oo. The disconnected 
contribution for the taste-axial- vector meson is 



and there is a similar contribution for the taste- vector meson. 

It is convenient to carry out a partial fraction expansion of the disconnected contributions 
as follows: 

rh^ ( h\rh^ (h\\ - ^r,sSr',s'SfJr,g' ( 3/2 1/2 \ 

KA-k)ct^fr',g'Ak))^^^ - — l^p^Mi; " PTMjj ^^^^ 

and 



{Kjri-k)<l>f'r'Mk)),,,, = -Sr,s5r',s'5f,5r.>5A + ^^^^ + p + ^2 



(49) 
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where 



Ml A - Ml A 


{M^A 




MhA) 




MIa - M'^A 




{M'uA 




M',a) 




MIa - M% 




(M'uA 




MI,a) 



(50) 



In the language of Refs. [13, 14], gu, g^,, and gr,' are simply the residues for Eq. (47). Similarly, 
the factors of 3/2 and —1/2 in Eq. (48) are the residues for Eq. (46). 



C. Isovector uq correlator 



We now specialize to the isovector oq correlator. We consider, for simplicity, the ud flavor 
state. Only the quark-line-connected contribution appears in the QCD correlator 



Baoip, ra) = Bu^Mp, ra) = - J^i'T exp(zp ■ x) (Tt[M-\x, r; 0, 0)M^'\x, r; 0, 0)]) . 



In terms of the meson fields, the bubble correlator is (for ra = tA) 



Bu^-uAiiP^i^) 



y r,s,f,br',s',f',b' 

4'urjs 



After carrying out the Wick contractions and switching to momentum space we get 



Baoip) 



A2(L/2)3(iV,/2) 

1 



- 4n,5] 

k 

- 4n^^ 

k 

- 4n^^ 



f,b k 
1 



+ M^^Ak + P? + M%^, 



SI 



[k + p)2 + Ml, J 3nr (P + M2,)(P + M^j) 

46a k^ + MIa 

{k + pf + Mi^ (F + Mi^W + M,2^)(P + M2^) 
45y e + Mly 



(51) 



(52) 



(53) 



(A; +p)2 + Mljy (F + M2^)(F + M2^)(P + 

Notice, in particular, the negative weight threshold in the second term and the spurious taste 
nonsinglet vrr/ thresholds involving Goldstone-boson-like members of the rj taste multiplet. 
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In the continuum limit, in which taste-symmetry is restored, we have 



BaM 



lA\L/2f{N,/2) 1 



k 



k^ + M%{k + pY + M% 



_(A; + p)2 + M2 (F + M2)(P + M2)J j ■ 
Here the total contribution from pairs of light states with mass Mu is proportional to 

(32n2 - 2) , 



(54) 



which vanishes when n,. = 1/4. The negative norm threshold has neatly canceled the 
unphysical thresholds. The surviving thresholds are the physical KK and taste singlet irrj. 



D. Isosinglet /o correlator 

In this case we use the isosinglet operator {puu,i + Pdd,i)/V^- We have both quark-line- 
connected and quark-line-disconnected contributions 



Bf^,{p,Ta) = 5/o,conn(P,™) + 5/0, disc (P,™) (55) 

Y^i-r exp(zp ■ xa) {Tt[M-\x, r; 0, 0)M-'\x, r; 0, 0)]) (56) 



-B/o,conn(P, ™) = -[Bu,u-,u,u,conn{P,T'a) + Bd,d-Ad,conn{p,Ta)] 



%o,disc(P,™) = -[Bu,u;u,u,disc{P,Ta) + B 

u,u;d,d,disc 

{p, to) 

+ Bd,d-u,u,disc{P, to) + Bd,d;d,d,disc{p, ™)] 

= ^Y.^Mw-xa){TT[M~\x,T-x,T)]TT[M-'\0,0-0,0)]) . (57) 

X 

The weight of the disconnected part is n//4 for Uf = 2 degenerate flavors for the state. The 
connected part of the correlator is identical to the full Cq correlator. 
In terms of the meson fields, the bubble correlator is (for ra = tA) 



Bfo{p,tA) = ^^exp(^p-yA) E 

y r,s,f,br',s',f',b' 



A^ 



4'ur,fs 

+ (0d.,/s(^)0?.A(^)'^l'J'«'(O)4«',*''(O)>] (5^ 
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In momentum space the correlator becomes 

Bfoi'P) 



A2(L/2)3(iVi/2) 
1 



f,b k 



b k 



- ^nrJ2 

k 

- 4n,^ 

k 

- 4n,^ 

k 

+ ^nlY. 

k 



1 



1 1 



(59) 



e + Mil 



45, 



A;2 + Ml,, 



[i^k + pY + (P + Mly){k^ + M^^y){k^ + M2^) 
1 {k + pf + Mlj 1 



M 



SI 



Mr [{k + p)2 + Mljj\[{k + p)2 + (P + M2,)(P + M2,)_ 

45aP + p)2 + M|^] 



X 



P + pY + M2^] [(A: + pf + M2^] [(A: + pf 
SA{k' + Ml^) 



(P + M2J(P + M2J(P + M2 )J 



k 



A5v[{k + pf + Mly] 



X 



P + p)2 + M2 ^] [(A: + p)2 + M2^] [{k + p)2 



^2 . M|^^ 



In terms of valence quark world lines the first five terms are quark-line connected and the 
last three are disconnected. 



In the continuum limit we have 
BfM 



A2(L/2)3(iV,/2) 
1 



+ 32^2^ 

k 

- 4n,^ 



16-.^ EE 

/ k 
1 



k^ + M%{k+pY + M% 



k'^ + M^ik + p^ + M^ 
1 1 



(60) 



P + M2 



+ Ant 



[k+pY + Ml 2,nr (P + M2)(P + M2) 
{k + pf + Ml 1 



k^ + Ml 



■in,, [{k + pf + Ml][{k + pY + M2] 3n,, (P + M2)(P + M2)_ 
The two-pion threshold {p + kY + M^ = and k^ + = has a weight proportional to 

(64n2 - l)/i2. (61) 

When 77,^ = 1/4 the weight is 3 (for three physical pion channels). Thus, once again, only 
physical thresholds survive the continuum limit. 
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E. Single-flavor staggered fermions 



Single-flavor QCD has no Goldstone bosons. The low-lying pseudoscalar (call it the r]') is 
lifted by the anomaly. With the staggered fermion action, however, only the taste-singlet rj' 
is lifted by the anomaly. The other 15 members of the taste multiplet (call them 77) remain 
light. The member with pseudoscalar taste is an exact Goldstone boson. Such a spectrum 
would seem to spell trouble for the rooted theory. It is interesting to examine the scalar 
meson (/o) correlator to see how the corresponding rooted chiral theory heals itself in the 
continuum limit. 

Call the single replicated flavor u. The connected meson correlator is as before [Eq. (45)]. 
We choose not to decouple the taste-singlet 77' in this case because it is the only physical 
meson. The disconnected correlator for the taste singlet is then 



f',g' 



disc 



1 



(62) 



Similarly, the disconnected correlators for the taste axial vector and taste vector can be 
written as 

-Sa 



disc ^r,sSr',s'SfJf',g' ^ M^^){k^ + M^, ^] 



(63) 



and {A ^ V), where in this case M^/^ = M^^ + rirSA, and similarly for M^, 



With these changes the /o correlator becomes 



1^' 



1 



1 



F + M2JA; + p)2 + M2, 
1 1 



+ 4n^E 

k 
k 

+ 2n^E 

k 

+ 2n^E 

k 

+ 2n^E 



[{k + py + M2^](P + M2^)(P + M^^) 
Sv 



(64) 



k 



{k + py + Ml J {k + pf + M2 , J Hr + Mil F + M2 ^_ 



[{k+pY + MlMk+P? + + Ml^W + M,%) 

Sy 

[[{k + pY + Mly][{k+pY + M2^](P + M2^)(F + M^^,y) 
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We note that a simplified version of our result (setting the discretization corrections from 
5a and 5y to zero) was presented previously [2]. In the continuum limit the would-be- 
Goldstone thresholds become degenerate with the negative-norm threshold, with a net weight 
proportional to 

(32^2 - 2) . (65) 

When rtr = 1/4, the would-be Goldstone bosons decouple from the /o correlator, leaving 
only the physical high-lying rj'T]' channel. 

IV. SIMULATIONS AND RESULTS 

In this work we analyzed the 0.12 fm ensemble of 510 24^ x 64 gauge configurations 
generated in the presence of 2 + 1 flavors of Asqtad improved staggered quarks with bare 
quark masses arriud = 0.005 and arris = 0.05 and bare gauge coupling 10/ = 6.76 [17]. 

We set valence quark masses equal to the sea quark masses. Table I gives the pseudoscalar 
masses used in our fits with the exception of the masses rjA, v'a^ Vvy Vv- Those masses vary 
with the fit parameters 6a and 6v- 

For the light quark Dirac operator M„, we measured the point-to-point quark-line con- 
nected correlator 

Cconnip, r) = Y^i-r cos(p ■ x) {Tt[M-\x, t; 0, 0)M-'\x, T- 0, 0)]) (66) 
and point-to-point quark-line disconnected correlator 

Cdisc(p, r) = Y.{-r cos(p ■ x) (TrM-i(x, r; f , r)TrM-i(0, 0; 0, 0)) . (67) 

X 

In the latter case we use noisy estimators based on random Z{2) color vectors [26] r]k for 
k = l,...iV = 200: 

TrM-i(f,r;f,r)TrM-i(0,0;0,0) ^ J_ E nk{S,T)M~\x,T-y)r^k{y) (68) 

^) k^k',y,y' 

X r]kiO,0)M-\0,0;y')rjk'{y'). 
In terms of these correlators the Oq and /o correlators are 

CfoiP.r) = Cconn{p,T) - ^Cdisc{p,r) (69) 
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Correlators in each channel were measured at five momenta p = (0, 0, 0), (1, 0, 0), (1, 1, 0), 
(1, 1, 1), and (2, 0, 0). All ten correlators were then fit to the following model 

Cao (P, t) = Cmcson,ao (P, t) + Ba^ (f, t) 

Cfo {P, r) = CmesonJo (P. t) + ^/o (P. t) (70) 

where 

Cmeson,ao (P, t) = Ko{p) GXp [--Ea^ (p)r] + K^A{p){-y CXp (p)t"] + (t ^ A^j - t) (71) 

Cmeson,/o(P7 ^) = Co(p) + &/o(p) exp[-E/o (p)r] + ^,a(p)(-)^ cxp (p)^] + (t A'f - r). 

This fitting model adds explicit and /o poles, as well as the corresponding negative parity 
states, to the bubble contribution. Such states are outside the scope of the low order chiral 
Lagrangian in Eq. (2). Of course it is possible to enlarge the Lagrangian to include them 
[21]. Taste-breaking effects complicate this exercise. Moreover, we would need to introduce 
a variety of higher order chiral couplings, which are unlikely to be well constrained by our 
data. Therefore, we took the more modest approach and treated these additional terms 
empirically, keeping in mind the possibility of higher order chiral effects. 

Our parameterization of the momentum dependence of the overlap factors hj{p) requires 
some discussion. The and /o are produced through the scalar density with spin-taste 
assignment 1x1. Thus at zeroth order in the ag — tt — r/ coupling their contributions should 
be inversely proportional to their energies 6j(0) = l/2Ej{p). At higher order an iteration 
of the bubble contribution alters the momentum dependence of the pole residue [21]. For 
present purposes we chose the empirical fitting form 

bj{p) = b,o + b,ip^. (72) 

and adjusted the constants bj^ and bji. 

The negative parity states are the taste-axial-vector pion tt^ and the taste-axial-vector 
7] A- As staggered partners to the ao and /o they couple through axial vector currents with 
spin-taste assignment 7075 x 7075, which contribute a factor of the energy to source and 
sink. Thus their bare momentum dependence should be proportional to their energies 

b,{p) = bjE{p). (73) 

We kept this form, adjusting bj. 
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The constant co{p) is zero for all momenta except p = (0,0,0), in which case it gives 
the vacuum-disconnected part of the /o correlator. There are eleven fit parameters for the 
meson terms alone, but the two negative parity masses were constrained tightly by priors: 
the tta, to the previously measured value, and the rj^, to the same derived mass that we 
used in the bubble term. 

The bubble terms Bag and Bf^ in the fitting function Eq (70) are given in momentum 
space by Eqs. (53) and (60). Their time- Fourier transforms yield Bagip, r) and Bf^^p, r) by 
applying the following identity term by term: 

1 p-ipot p-[E^{k)+E2(k)]t 

Biv t) oc = V - = V - f74) 

^ A\N,/2f^,{k^ + Ml)\{p-kf + M',] ^ AE,{k)E,{k) 



where Ej{k) = y + Mj, and, as usual, tA = ra. Thus, for example, the KK contribution 

to Bag {p, t) for taste h is 

2 -[EKb<,k)+EKb[k)\t 

The bubble terms Bag{p,T) and Bf^^{p,T) were parameterized by the three low energy 
couplings /i = m^/(2m£), 6a = 0^5^, and 6v = a^S'j^^ in the notation of Ref. [15]. They were 
allowed to vary to give the best fit. The taste multiplet masses in the bubble terms were 
fixed as noted above. The sum over intermediate momenta was cut off when the total energy 
of the two-body state exceeded 1.8/a or any momentum component exceeded 7r/(3a). We 
determined that such a cut off gave acceptable accuracy for r > 4. 

In summary, we fit all ten correlators with fourteen parameters, eleven of which were 
needed to parameterize the four explicit meson terms and three low energy couplings were 
needed for the bubble contribution. Through a prior, we constrained the value of 6v to 
conform to previous fits to the pseudoscalar masses and decay constants [15], leaving only 
two of the low energy couplings to be adjusted independently. Our best fit gave x^/dof = 
126/109 (CL 0.13). 

The fitted functional form is compared with the data in Figs 1-3. 

Results of the fits are compared with results from fits to the meson masses and decay 
constants in Table II. The agreement is worse if we used the bare value rifi = 4.5 from 
those fits, rather than the higher-order m^/(2m^), suggesting, perhaps, that a higher order 
calculation of the bubble contribution might improve the agreement. 
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FIG. 1: Best fit to the oq correlator for five total cm momenta. The fitting range is indicated by 
points and fitted lines in red and blue (darker points and lines). Occasional points with negative 
central values are not plotted. 

Our fit Meson masses and decays 
rim2/(2m,,d) 7.3(1.6) 6.7 
6v (prior) -0.016(23) 

6a -0.056(10) -0.040(6) 

TABLE II: Comparison of our fit parameters for the rS^PT low energy constants with results from 
[15] 

The fitted masses of the ao and /o in units of the lattice spacing are 0.61(5) and 0.45(9), 
respectively. 

V. SUMMARY AND CONCLUSIONS 

We have derived the two-pseudoscalar-meson "bubble" contribution to the /q correlator 
in lowest order S^PT, thereby extending the result for the Cq in reference [20]. We have 
used this model to fit simulation data for the point-to-point and /o correlators and found 
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FIG. 2: Best fit to the /o correlator for four total cm momenta. 
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FIG. 3: Best fit to the zero momentum /o correlator. 

that best-fit values of the three chiral low energy couplings are in reasonable agreement with 
values previously obtained in fits to the light meson spectra and decay constants [15]. 
The two-meson bubble term in S^PT provides a useful illustration of the lattice artifacts 
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induced by the fourth-root approximation, since it involves quark loops coming from the 
fermion determinant. The artifacts include thresholds at unphysical energies and thresh- 
olds with negative weights. These are the same sorts of artifacts commonly observed with 
quenching or partial quenching. These contributions are clearly present in the Oq and /o 
channels in our QCD simulation with the Asqtad action at a = 0.12 fm. We have found 
that they must be taken into account in a successful spectral analysis. Fortunately, rS^PT 
provides an explicit parameterization of their contributions for the interpolating operators 
we have chosen, thereby allowing a fit to simulation data with a manageable number of 
parameters. The rS^PT predicts further that these lattice artifacts disappear in the contin- 
uum limit, leaving only physical two-body thresholds. This result is in full accordance with 
the fourth-root analysis of Ref [2]. It will be interesting to see whether this expectation is 
borne out in numerical QCD simulations at smaller lattice spacing. 
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